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Abstract 

We explore a possibility for the Universe to cross the w = —1 cosmological constant 
barrier for the dark energy state parameter. We consider the Universe as a slowly 
decaying D3-brane. The D3-brane dynamics is approximately described by a nonlocal 
string tachyon interaction and a back reaction of gravity is incorporated in the closed 
string tachyon dynamics. In a local effective approximation this model contains one 
phantom component and one usual field with a simple polynomial interaction. To 
understand cosmological properties of this system we study toy models with the same 
scalar fields but with modified interactions. These modifications admit polynomial 
superpotentials. We find restrictions on these interactions under which it is possible to 
reach w = —1 from below at large time. Explicit solutions with the dark energy state 
parameter crossing/non-crossing the barrier w = —1 at large time are presented. 
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1 Introduction 



The combined analysis of the type la supernovae, galaxy clusters measurements and WMAP 
data provides an evidence for the accelerated cosmic expansion PQ 121 El • The cosmological 
acceleration strongly indicates that the present day Universe is dominated by smoothly 
distributed slowly varying Dark Energy (DE) component. The modern constraints on the 
DE state parameter are around the cosmological constant value, iy = — liO.ljHll^ and a 
possibility that w is varied in time is not excluded. From the theoretical point of view there 
are three essentially different cases: w > —1 (quintessence), w = —1 (cosmological constant) 
and w < —1 (phantom) (0101111 and refs. therein). 

Since from the observational point of view there is no barrier between these three pos- 
sibilities it is worth to consider models where these three cases are realized. Under general 
assumptions it is proved in [8J that within one scalar field model one can realize only one 
possibility: w > —1 (usual model), or w < — 1 (phantom model). It is interesting that the 
interaction with the cold dark matter does not change the situation and does not remove the 
cosmological constant barrier jH]. There are several phenomenological models describing the 
crossing of the cosmological constant barrier JUj- Most of them use more then one scalar 
field or use a non-minimal coupling with the gravity, or modified gravity, in particular via 
the brane-world scenarios. In two-field models one of these two fields is a phantom, other 
one is a usual field and the interaction is nonpolynomial in general. 

It is important to find a model which follows from the fundamental principles and de- 
scribes a crossing of the w = —1 barrier. 

In this paper we show that such a model may appear within a brane approach when the 
Universe is considered as a slowly decaying D3-brane and a possibility to cross the barrier 
comes from taking into account a back reaction of the D3-brane. This DE model jHI assumes 
that our Universe is a slowly decaying D3-brane and its dynamics is described by the open 
string tachyon mode and the back reaction of this brane is incorporated in the dynamics 
of the closed string tachyon. The open string tachyon dynamics is described within a level 
truncated open string field theory (OSFT). The notable feature of this OSFT description 
of the tachyon dynamics is a non-local polynomial interaction ^TJ-PH]- It turns out the 
open string tachyon behavior is effectively described by a scalar field with a negative kinetic 
term f phantom) |16| IT7j . However this model does not suffer from quantum instability, which 
usually phantom models have, since in the nonlocal theory obtained from OSFT there are 
no ghosts at all near the non-perturbative vacuum j^l. 

Level truncated cubic OSFT fixes the form of the interaction of local fields to be a cubic 
polynomial with non-local form-factors. Integrating out low lying auxiliary fields one gets a 
4-th order polynomial Higher order auxiliary fields may change the coefficients in front 
of lower terms and produce higher order polynomials. All these corrections are of higher 
orders of a' . 

The second scalar field comes from the closed string sector, similar to and its effective 
local description is given by an ordinary kinetic term ^H] and, generally speaking, a non- 
polynomial self-interaction [201 ■ An exact form of the open-closed tachyons interaction is 
not known and we consider the simplest polynomial interaction. 

Our goal is to understand the following: is it possible in the two component polynomial 
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model that w crosses the barrier w = —1 at large time and reaches to —1 from below at 
infinity? For this purpose we study special polynomial two component models. For these 
models there exist third order odd superpotentials. An existence of a superpotential puts 
restrictions on the form of the potential. For polynomial potentials these restrictions give 
relations among coefficients. In this polynomial case we can estimate the behavior of DE 
state parameter at the large times. We expect that small variations of the coefficients of the 
potentials obtained from the given superpotential do not change qualitatively a behavior of 
the system. 

The superpotentials under consideration produce the potentials which are rather close 
to the form of the open-closed tachyon potential for a non-BPS brane. Indeed, within the 
level truncated sting field theory description of a non-BPS D3-brane decay both fields have 
tachyon mass terms and the interaction is polynomial, the 4-th order at the lowest levels. 
A natural deformation of this form of the open-closed string tachyons potential is given by 
extra 6-th order terms. 

Corresponding local models in the fiat background admit exact solutions. An exact solu- 
tion of effective local model describing the pure open sector of a non-BPS D3-brane is given 
by the kink solution [T^] and the closed tachyon dynamics under reasonable assumptions is 
given by a lump solution 1^ • In a non-flat background there is a deformation of the 
effective local model describing the pure open sector of a non-BPS D3-brane such that the 
corresponding Friedmann equations have exact solutions [7]. A more straightforward gener- 
alization of the model [Zj to the case of two flelds gives a model with a kink-lump solution. 
This solution at late times has a behavior as a quintessence model, i.e. w goes to —1 from 
the above. 

We also construct an exactly solvable stringy DE model with the state parameter which 
crosses the cosmological constant barrier k; = — 1 at a rather late time from the above, reach 
its local minimal values that is less then —1 and approaches —1 from the below at inflnite 
time. The form of the potential in this case is rather complicated and we cannot construct 
it from the string fleld theory yet. The Hubble parameter in this model is a non-monotonic 
function of time as well as the DE state parameter w. 

2 The Model 

We consider a model of Einstein gravity interacting with a single phantom scalar fleld (j) and 
one standard scalar fleld ^ in the spatially flat Friedmann Universe. Since these scalar flelds 
are assumed to come from the string fleld theory the string mass Mg and a dimensionless 
open string coupling constant go emerges. In the typical cases phantom represents the open 
string tachyon and the usual scalar fleld the closed string tachyon jHl El 1^ ■ The action is 

S = I d'x^ (^-^R+ (^+^g^''d,<j>d,<j> -^g^'^d,^^^^ (2.1) 

where Mp is the reduced Planck mass, gfj_i, is a spatially flat Friedmann metric 

ds^ = - dt^ + o?{t){dxl + dxl + dxl). 
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and coordinates {t, Xi) and fields and ^ are dimensionless. Hereafter we use tfie dimension- 
less parameter rrip for short: 

m^ = ^ (22) 

s 

If the scalar fields depend only on time then equations of motion are as follows 

2H = \U''-iA, (2.3b) 
V / 

dV 

(P + 3H(P = — , (2.3c) 
dV 

^ + 3H^ = -—. (2.3d) 

Here dot denotes the time derivative and H = a(t)/a(t). 

The form of the potential is assumed to be given from the string field theory within the 
level truncation scheme. Usually for a finite order truncation the potential is a polynomial 
and its particular form depends on the string type. 

In the present analysis we impose the following restriction on the potential: 

• potential admits an existence of a polynomial superpotential (see details in and in 
the next section) 

• potential is even 

• 0(t) has non-zero asymptotics and ^(t) has zero asymptotics as t ^ oo 

• potential is not more than 6-th power 

• coefficient in front of 5-th and 6-th powers are of order l/m^ and the limit oo 
gives a nontrivial 4-th order potential. 

Particular exact solutions will be found by using more specific anzatses. We will see 
that for the solution to be constructed in Section E] the form of the potentials in the limit 
rrip — s> oo reproduces the one given by an approximation of the lowest level truncated string 
field theory. 



3 w = —1 Barrier for Two-component model with Poly- 
nomial Superpotential 

3.1 Setup 

We can assume that H{t) is a function (named as superpotential, see for example |22]) of 
(j){t) and ^{t): 

Hit) = Wi^it),m)- 
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This allows us to rewrite ()2.3b|) as 



dW ■ dW ■ I f 12 /O.N 



The latter system is certainly solved provided the relations 

dW 1 • 



d(f) 

dW 1 



^ (3.2) 



2m^ 



are satisfied. If this is the case we have the following relation between the potential V and 
the superpotential W 

This relation gives the potential in terms of W and its first derivatives with respect to and 
^. Provided the superpotential is given to find a solution of the dynamical system one has 
to solve the second order system of partial differential equations (|3.2|) . 

3.2 Construction of the potential 

In this subsection we construct polynomial potentials that admit a polynomial superpoten- 
tial. Recall, that we restrict ourself to have the maximal power in the potential to be equal 
to 6 and the potential should be even. Then general substitutions for <p{t) and ^(t) are as 
follows 

m,n=0,l,2 

Equivalence of second mixed derivatives of W implies 

2^12 = -P21, Xu = -2p20, a^Ol = -PlO, Pll = -2X02, X21 = Pl2 = X22 = P22 = 0. 

For the potential to be even we have to set to zero constants poi, pio, xqi, xiq and an 
integration constant in the W should be zero as well. Also in order to have the maximal 
power 6 in the interaction potential for (p and ^ we have to put p2i = and Xi2 = 0. 
Substituting expressions (j3.4j) into (j3.2|) after integration we have 

W = (^Poo4> + ^P02(P^ - Xoo^ - ^^20^^ + P2O0^^ " Xo2(p^^^ ■ (3.5) 
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One can obtain the form of the potential V from the relation (j3.3|) . However, we postpone 
this to the next subsection when asymptotic late time behavior will be specified. 

Note that in the case of one field the superpotential W defines this scalar field, as a so- 
lution of the first order differential equation, which always can be trivially solved in quadra- 
tures |22] and there is no difference to start with explicit form of the (phantom) scalar field 
as a function of time or the corresponding form of superpotential. In the case of two fields 
the superpotential method gives the second order system of differential equations, which may 
be non-integrable. In this case it is more preferable to start from the form of superpotential, 
which corresponds to the required form of the potential. In Section El we demonstrate that 
the scalar and phantom scalar fields can have very unusual dependence on time, which can 
not be predicted from the consideration of models with one field and a polynomial potential. 



3.3 Time evolution 

Differential equations for our fields when all relations among p^n and Xmn constants are 
taken into account read as follows 

= Poo +PO20'(t) - 2Xo20(t)eW +P20e'(t), 

To specify the boundary conditions let us recall that we have in mind the following picture. 
We assume that the phantom field smoothly rolls from an unstable perturbative vacuum 
{(f) = 0) to a nonperturbative one, say, = a and stops there. The field ^ we expect to go 
asymptotically to zero in the infinite future. An asymptotic behavior implies poo = — Po20^ 
and Xoo = —xo2Ci^ and we left with the following system 



m = -P02a' +PO20'(t) - 2X02<P{t)^{t)+p20^\t), 

= -xo2a^ + a;o20'(t) - 2p2o0(t)e(^) + X2o^\t). 
The superpotential W can be rewritten in the following form 



(3.7) 



^ = A (-Po20(3a2 _ ^2) ^ 3p^^^^2 ^ Sxo2aa' - 0') - X2oe) ■ (3-8) 
brrip 

The corresponding potential V is the following 

V = ^{-po2a^ + Po20^ - 2xo20^ + P2o^^f - ^{-xo2a'^ + xo2(j)^ - 2p2o0^ + 2:20^^)^+ 



+ (Po20(3a2 _ ^2) _ 3^Q2e(a' - 0') - 3^200^' + 0:20^')' 



(3.9) 



p 



3.4 Cosmological consequences: late time behaviour 

From the cosmological point of view we address the following questions to our model. What 
is the behavior of the Hubble parameter H, how the state parameter w and the deceleration 
parameter q do evolve? 
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Even without having a time dependence of the fields and ^ we can answer some of the 
above question provided that we know the asymptotic behavior of the fields. Indeed, we 
assume the field starts from and goes to a finite asymptotic a and its velocity goes to 
zero in the infinite future. The field ^(t) and its velocity ^{t) go to zero in the infinite future. 
Recall, that the function H{t) is restored once we substitute the time dependence and 
^(t) into ()3.8|) . As the first result we see that H{t) in the infinite future goes asymptotically 
to the following value 

(^•^'') 

We immediately see that po2 should be negative if a is a positive asymptotic of the field 
Also it is evident that H{t) goes to zero. 

Further one can expand functions (f){t) and ^{t) for large times as follows 

0(t) = a + /(t) + ..., m = 9{t) + - (3.11) 

where f{t) <C a, g{t) <C a and the ratio f(t)/g{t) is finite. Assuming such an expansion we 
have the following asymptotic behavior of the Hubble parameter 

= + {po2f{t) - 2x,,f{t)g{t) + p,o9\t)) . (3.12) 

3mp imp 

The eigenvalues of the quadratic form in / and g are found to be 

A//,i = ^(P20 + Po2 + A/(P2o-Po2)^ + 4a;o2), 

Aj/,2 = ^(P20 + P02 - \/(P20 -P02)^ + 4X02) 

and they determine whether H{t) comes to its asymptotic value from the above {\h,i > 
and Xh,2 > 0) or from the below {Xh,i < and Xh,2 < 0). If Xh,i and Xh,2 have opposite 
signs we need to use more detailed approximation. 

Now we turn to the behavior of the state parameter and the deceleration parameter. 
They are related with the Hubble parameter by the following relations 

2 H H 
^ = -1-3^' ^ = -1-^- 

Since H{t) in our consideration goes asymptotically to a finite constant and its time derivative 
vanishes both the state and deceleration parameters go to —1. The question does w approach 
to —1 from the above or from the below is very important. The first case is the so called 
quintessence like behavior and the second case is the phantom like behavior. It is convenient 
to rewrite the relation for the state parameter using the equation ()2.3b|l as follows 
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Substituting expressions for the and ^(t) from (j3.7p we get 



where 



(3.14) 



We employ again the asymptotic expansion ()3.11|) to write 

Aas = ^a\pl^ - xl^)f{t) - Aa\plo - 42)9^^) + 80^X02(^20 - Pm)fit)g{t). (3.15) 
The quadratic form 

(Po2 - fit) - (pIo - 42)9\t) + 2a;o2(p2o - Po2)fit)git) (3.16) 
has the following eigenvalues 

^;.^,-.^. + ./(P^, + P^.p + 4.,,(.,,-....o2), 

(3.17) 



Aa,i = ^(Po2 -plo + y (Po2 + pIoY + 4^^2(42 - 2P20P02), 



Aa,2 = -(Po2 - P20 - V (P02 + P2o)^ + 4X^2(2:^2 - 2P20P02)- 

Therefore, when both Aa,i and Aa,2 are positive we have a phantom like behavior, when 
these A-s are both negative we have a quintessence like behavior. When Aa,i and Aa,2 have 
opposite signs we may have oscillations at large times near the cosmological constant barrier 
w = —1. 

In the next sections we consider two special solutions. The first one corresponds to the 
quintessence behavior and the second one to the phantom behavior. Moreover we will see 
that for these solutions the state parameter crosses the w = —1 barrier. Notice that such a 
crossing is forbidden in one field models 

4 Quintessence late time solution 
4.1 Anzats and corresponding potential 

We are about to construct a solution to the system (j3.7j) . The system is essentially simplified 
if we take 

X02 = X20 = 0. (4.1) 

The latter is the anzats we explore in this section The solution which will be found possesses 
the properties reflecting its name. Substitution of this anzats into ()3.7|) gives 

0(t) = -po2a^ + Pm<f (t) + P2o^^ (t) , 

m = -2p2,mm- 



The superpotential W given by (j3.5|) for the case (j4. 1|) reads as follows 
The corresponding potential V can be found using the relation ()3.3|) to be 



V = l {-P02a^ + Po20' + P2o^^Y - 2plo4>^^^+ 



2 

+ ^0' (Po2(3a2-02)_3p2^^2)^ 



(4.4) 



4.2 Solution 



A solution to the system ()4.2j) when filed starts from and goes asymptotically to a and 
field ^ asymptotically vanishes is the following 

= atanh(2ap2o^) (4-5) 

and 

e = . (4.6) 

cosh(2ap2o^) 

Hereafter in this section we denote r = ^02/^20- 

Let us note that one obtains the same solution ()4.5|1 . ()4.6|) for different potentials. Namely, 
the solution is not violated if we take a potential of the form 

Vi = V + 6V, (4.7) 

where V is the potential given by ()4.4|1 and 6V is such that SV, d{6V)/d(f) and d{6V)/d^ 
are zero on the solution. For 0(t) and ^(t) given by ()4.5j) and ()4.(i|l respectively the most 
general form of even 6V with the maximal power of interaction 6 is the following 



6V = A 



1 ' 
2 + r 



This example shows that the same functions (p{t), ^(t) (and consequently the Hubble 
parameter H{t), state parameter w and deceleration parameter q{t)) can correspond to 
different potentials V{(f),C,). 

4.3 Cosmological properties 

We obtain H(t) by substituting ()4.5|) and ()4.6|) into ()4.3|) and expressing the result as a 
function -/^(t). The result is 
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The function H{t) is not monotonic for general values of the parameters and has an extremum 
at the point 

tc = ^4 • ^-lO 

2ap2o 

We are certainly interested in the case when this tc is real, i.e. the argument of the 
logarithm should be a positive real value. That means that we have to have r > —2. 
Moreover, if r > —2 then the argument of the logarithm is greater than 1, and consequently 
the value of the logarithm is positive. Further we recall that in order to have a positive 
asymptotic for H{t) we have required po2 < 0. On the other hand expression for the tc implies 
that P20 > if we are interesting in positive time semi-axis (the situation is symmetric for 
the negative time semi-axis). Thus, r turns out to be less than 0. Eventually, we state that 

-2 < r < 0. 

The Hubble parameter in the extremum is 



2a^ I 1 



Recall that at large times the Hubble constant goes to 



a3 



Hoc = -T. 

omi 



and the ratio Hd is as follows 



-f^oo r V 3 + r 



^ I ' (4.12) 



and is determined by the ratio r of parameters po2 and p2o- It is a matter of a simple algebra 
to check that for r > —2 the ratio (j4.12p is greater than 1. This means that for a specified 
domain of r the point tc corresponds to a maximum. The typical plots corresponding to the 
performed analysis are shown in Fig. ^ In the case r < —2 the tc becomes imaginary and 
the function H{t) turns out to be monotonic. This situation is close to the one field model. 
The case r > is implausible because H{t) changes the sign during the evolution and has a 
negative asymptotic. 

The state parameter w is given by the following expression 



w = —14 



X 



a2 cosh^(2ap2oi) - P02 - 

1 — 3 tanh^(2ap2oi) 2po2 



2,. ,x !! X (4-13) 



cosh {2ap2ot) - 1 P02 cosh {2ap2ot) - po2 - 3p; 



'20 



It has a singularity in the origin and behaves as —1/t^. At the point tc the state parameter 
w crosses —1 because at this point H{t) = 0. After tc for particular values of the parameters 
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Figure 1: The Hubble parameter H{t), total w and total q for = 0.2, a = 1, p2o = 1/2 
and po2 = —0.01 

(see Fig. [T]) appears a period of deceleration (g > 0), however, at late times the Universe 
returns to the acceleration, w and q for this solution approach —1 from the above. The latter 
is evident from the expression for A ()3.14|) 



A 



-4^20 a 



2 + r 



3 + r 



cosh (2ap2o^) cosh (2ap2o^) 



For the large time the first term in the parentheses dominates and since in our case it is 
assumed that r > — 2 we obtain that A < 0, w goes to —1 from the above and the solutions 
have the quintessence like behavior. 



4.4 Connection to SSFT 

The potential (j4.4|) contains mass terms for the fields 
follows 



and ^. Their masses are given as 



2 2 



nil: 



3a^ 
2m^ 



-2po2P2oa 



(4.14) 



However we have obtained previously the following restrictions: po2 should be negative, p2o 
should be positive once the asymptotic a is chosen to be positive in order to have a suitable 
cosmological behavior. Also it follows from ()4.12j) that the ratio r should be small if we 
want to observe a large ratio of the maximal Hubble constant and the asymptotic Hubble 
constant. These restrictions say us that both and are small, the field is a tachyon 
in the limit of the large reduced Planck mass and the field ^ has a positive mass squared. 

The situation drastically changes once we make use of the freedom ()4.8j) . We can choose 
for simplicity f i = f2 = fs = 0. In this case new potential will give new masses for the fields. 
Indeed, 



Mi 



2 2 
P02a 



Ml 



3a^ 
2m^ 



-2po2P2oa - 



2 - 4Aa2, 



AAa^ 



(4.15) 



2 + r' 
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Provided po2 is small we effectively change the character of the fields because now if A > 
they are both tachyons. Moreover, in the limit r —>■ the ratio of masses M^/M^ goes to 2 
as it should be if we consider as the open string tachyon and ^ as the closed string tachyon. 
The trajectories of fields and ^ are presented in Fig. |2l (red lines). 




Figure 2: Potential ()4.4j) for = 0.2, a = 1, p2o = 1/2 and po2 = —0.01 (left) and potential 
(j4.7|) for p = a = T = 0, Aa =10 and all other parameters the same as in the left plot 
(right). 



5 Phantom late time solution 
5.1 Anzats and corresponding potential 

In this section we construct a more involved solution to the system (j3.7|) which exhibits a 
phantom like late time behavior for particular values of the parameters. Let us assume the 
following relations among the coefficients of the system ()3.7p 

P02 + X02 = -ixo2 + P20) = P2Q + X20 = "C. (5.1) 

Under this assumption we have 

^{t) = -po2a' + Po2^\t) + 2(po2 + c)4>{mt) + (P02 + 2c)e{t), 

m = {Po2 + c)a^ - {Po2 + c)0'(t) - 2{po2 + 2c)0(t)e(t) - {po2 + 3c)e'(t). 

The superpotential under the conditions ()5.1|) has the form 
W=-^ (-Po20(3a2 _ ^2) ^ 3^^^^ ^ 2c)0e' - 3(po2 + c)e(a' - </>') + {po2 + Sc)^) (5.3) 
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and corresponding potential is 

V = i(-po2a' + + 2(po2 + c)0e + (P02 + 2c)e')'- 



- 2^(^02 + c)a^ - {Po2 + c)(l)^ - 2{po2 + 2c)0e - {Po2 + 3c)e?+ (5.4) 
1 



+ Y^(^^o20(3a2 _ ^2) ^ 3(^^^ ^ ^)^^^2 _ ^2) _ 3(p^^ ^ 2c)0e' - {po2 + 3c)e')'. 
5.2 Solution 

Comparing two lines of the system ()5.2j) one readily finds 

m+i{t)=ca'-c{<P{t)+mY- 

Therefore 

+ ^(t) = atanh(ac(t + to))- 
Substituting C,(t) = atanh(ac(t + to)) — 4>(t) iiito (j5.2|) one finds 

a^(Po2 + c)t-C 
cosh^ (ac(t + to)) 

and 



a tanh(ac(t + to)) f 2! " . , (5-5) 



e = _'.Z (5-6) 



Q^(P02 + c)t - C 
cosh^ (ac(t + to)) 

A behavior of the solution depends on particular values of parameters a, c, to and C. We 
adjust C in such a way that 0(0) = 0. This gives 

C = — — sinh(2acto). 

The form of trajectories ()5.5|) and ()5.6|) for particular values of the parameters is presented 
in Fig. El 

5.3 Cosmological properties 

On solutions and ()5.(jp the Hubble parameter has the following form 

Hit) = tanh(ac(t + to)) (2 + ^ 

6m^ V coslr(ac(t + to)), 

2tac(po2 + c) + sinh(2acto) 



4m2 cosh (ac(t + to)) 



An analysis of this function is rather involved because we arrive to transcendent equations 
once we want to find the extrema. The situation is simplified a bit if we use the equation 
(I2.3b|) to express the H(t). One can write 

Hit) ~ 02(t) - i\t) = (0(t) - e(t))(0(t) + m). 
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Figure 3: (j){t) and ^{t) for a = 1, c = 0.55, po2 = -0.05, to = -2 and = 0.2. 



The last multiplier for our solution is equal to 



cosh {ac{t + to)) 



The latter expression has the same sign as c and becomes only in the infinite future. Note 
that c should be positive if a is taken to be positive. Otherwise the asymptotic value of H(t) 
will be negative. Thus, the zeros of the H{t) are determined by zeros of the first multiplier 
<J)[t) —C,{t). Also, the sign of the first multiplier can determine the late time behavior. Using 
exact dependence for the fields one can write 



0(t) - iit) = a- 



,2c(2a(po2 + c)t + sinh(2acto)) tanh(ac(t + to)) 



2p, 



02 



cosh (ac(t + to)) 



This expressions leads to the following consequences. The late time behavior is governed by 
the sign of the sum po2 + c. It should be positive if we expect the phantom like late time 
behavior. Also we observe that a natural choice to = does not lead to new interesting 
cosmological properties. Indeed, in this case the expression (j5.8p is governed by the function 
ttanh(act) which is monotonic at t > 0. Thus, the nominator of (j5.8|) has not more than one 
zero. If so, then we cannot have a large pick for the function H{t) and a phantom like late 
time behavior simultaneously (exactly such an evolution is interesting from the cosmological 
point of view) because to posses these properties H{t) should have a local minimum. To 
summarize we have to have a > 0, c > 0, po2 < 0, to should be finite if we want to observe 
new effects and po2 + c > for the phantom like late time behavior. 

The complete analysis of zeros of the equation ()5.8p is very cumbersome. However, it is 
possible to find particular values of the parameters for which it has two roots. This situation 
is demonstrated in Fig. (jH). 
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Figure 4: The Hubble parameter H{t) and its fine structure (top pictures). The state 
parameter w{t) and deceleration parameter q{t) (bottom pictures). Here a = 1, c = 0.55, 
po2 = -0.05, to = -2 and = 0.2. 

6 Conclusion and Discussion 

In this paper we have investigated the dynamics of two component DE models, with one 
phantom field and one usual field with special polynomial potentials. The main motivation 
for us was a model of the the Universe as a slowly decaying D3-brane which dynamics is 
described by a tachyon field jB] . To take into account the back reaction of gravity we take one 
more scalar field. This scalar field has usual kinetic term. The model is close to the model 
considered by j2S] which is also considered in the DE context. Note also that in the closed 
bosonic string sector an extra phantom appears Within two component DE models 

with a general class of interactions which correspond to polynomial superpotentials we have 
found restrictions that show whether the model is a phantom-like [w goes to -1 from below), 
or it is a quintessence-like {w goes to -1 from above). 

In particular, for the simplest model inspired by a D3-brane we have found that an 
inclusion of the closed string tachyon drastically changes the late time regime and for two- 
component model w > — 1 at large time, while in open string case w < —1 at large time. 

The model considered in this paper is close to the model considered in This model 
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is unstable, while a stability of our model is provided by its string origin. Note also that in 
the closed bosonic string sector appears an extra phantom [24j. 

The two-component model considered in this paper is interesting also by the following 
reason. There are several attempts to unify the early time inflation with late time accelerated 
universe (see for example and refs. therein). Generally speaking it is rather difficult to 
do this mainly because the ratio of the Hubble parameters in the end of inflation to its value 
during the period of the late acceleration should be very large. In our case we have such a 
possibility just by taking r to be close to in formula ()4.12|) . 

Let us recall that two scalar fields, both with usual kinetic terms, have been used in the 
hybrid infiation |26], and in particular in [2Zj the superpotential has a simple quadratic form. 

We have also found the superpotentials depending on two components for which w < —1 
at late times. We have presented the explicit solution realizing this possibility. It would 
be very interesting to study small deformations of the corresponding potential and to make 
clear does the constructed solution stable or not under deformations of the form of potentials 
and after including the CDM. 
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